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Abstract 

We show from first principles the emergence of classical Boltzmann equations for a generic quantum 
field that couples to a collection of background fields and sources. Our analysis is based on relativis- 
tic nonequilibrium quantum field theory as described by the Kadanoff-Baym equations and does 
not rely on the usual semi-classical approximations, i.e. on-shell condition, close-to-equilibrium 
assumption, gradient expansion and the Kadanoff-Baym ansatz. Instead, we find simple analytical 
solutions to the full Kadanoff-Baym equations by using the Wentzel-Kramers-Brillouin (WKB) 
method, and discuss under which circumstances Boltzmann behavior emerges. 
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I. INTRODUCTION 



Nonequilibrium phenomena play a crucial role in many areas of physics, including the 
early history of the universe, heavy ion collisions, condensed matter physics and quantum 
information. In the era of precision cosmology and with the arrival of the LHC and RHIC 
experiments, in particular the first two applications, which require a relativistic description, 
have gained considerable interest. Transport in nonequilibrium situations can often in very 
good approximation be described by Boltzmann equations (BEs). These assume that the 
system can be characterized by a number of distribution functions for classical particles, 
which propagate freely between isolated interactions and carry no memory of their history. 
However, the definition of asymptotic states, on which the single particle description is 
based, is ambiguous in a dense plasma. What is more, the standard BEs by construction 
cannot describe memory and off-shell effects or quantum coherence. Usually these effects 
are treated by effective kinetic equations of the Boltzmann type [l- 21] , i.e. by a set of first 
order differential equations for generalized distribution functions that are local in time. As 
the above issues are conceptual, their range of validity and possible corrections cannot be 
determined within a framework of BEs and require a derivation from first principles. 



The full equations of motion of nonequilibrium quantum field theory, on which first- 
principle derivations of the BEs usually are based, are known as Kadanoff-Baym equations 
(KBEs) [22]. These equations, being coupled second order integro-differential equations, are 
considerably more complicated than BEs. Analytical approaches to establish a connection 
between out-of-equilibrium quantum fields and kinetic equations usually make a number of 
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-]2l|): These include 



assumptions on the KBEs before they are solved (e.g. Refs. [lO, 
a gradient expansion and/or close-to-equilibrium assumption for all fields, the quasiparticle 
approximation and the Kadanoff-Baym ansatz. In this letter, we show how the full KBEs 
can be solved by using the WKB [231] method (for earlier uses of the WKB method in a 



similar context see e.g. Refs. [Hi, |24j). For the case of a single real scalar field, we find 
analytical expressions for its two-point correlation functions that are valid at all times, and 
discuss under which circumstances they exhibit Boltzmann behavior. The approach can be 
generalized to scenarios with multiple fields; backreaction can be included self-consistently. 



II. NONEQUILIBRIUM QUANTUM FIELD THEORY 



We consider the dynamics of a real scalar field that is described by relativistic quantum 
field theory. The field weakly couples to a background, possibly containing many degrees 
of freedom whose dynamics is in principle known and that we refer to as Xi- The Lagrangian 
reads 

c = ^oyo^ - \m{tf4> 2 - o[ Xi ] + C Xt , (1) 

where 0[x%\ denotes a generic combination of fields Xi-, an d jC Xi determines the dynamics 
of Xi ( we use h — c — 1). We allow a time- dependent mass m{t) to account for Hubble 
expansion when interpreting t as conformal time. 

In quantum physics, any thermodynamic system can be characterized by a density matrix 
g. Knowledge of the density matrix, or knowledge of all n-point functions (0(xi) • • ■ (p(x n )) 
of the fields, allows to compute expectation values for all observables at all times (with 
(•••)= tr[g •■■]). However, most quantities of practical interest for which one formulates a 
Boltzmann equation can be expressed in terms of one- and two-point functions; this includes 
the energy-momentum tensor and charge densities. It is, therefore, usually sufficient to track 
the time evolution of these. 

An out-of-equilibrium quantum field has two independent connected two-point functions. 
In case of they are conveniently chosen as 

A-(x 1 ,x 2 ) = i([<f ) (x 1 ),(t)(x 2 )}) , (2) 
A + (xi,a; 2 ) = i({0(^i),0(^ 2 )}) , 

with the obvious symmetry relations A ± (x 2 ,xi) = ±A ± (x±, x 2 ). Here [, ] and {, } are com- 
mutator and anti-commutator, respectively. A~(a;i,a;2) is known as spectral function and 
basically encodes information about the spectrum of resonances. The statistical propagator 
A + (xi,x 2 ) carries information about the occupation numbers of different modes and can be 
viewed as a generalization of the classical particle distribution function. We have in mind 
applications in cosmology and restrict the analysis to spatially homogeneous systems. Then, 
the correlation functions only depend on relative spatial coordinates Xi — x 2 etc., and it is 
convenient to perform a spatial Fourier transform in these coordinates, yielding functions 
like Aj(t!,t 2 ) = / d 3 ( Xl -^ 2 )e-^-^A ± (x 1 ,x 2 ). 
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In a general out-of-equilibrium system the two-point functions A^ (t±, t 2 ) have to be found 
as solutions to the KBEs 

(dl +u n (t 1 ) 2 ) A-(t!,t 2 ) = - f h dt' n-^^OA-^',^) (3) 

(^ 2 i +u; q (t 1 ) 2 )A+(t 1 ,t 2 ) = - t dt'Il-(t 1 ,t')A+(t',t 2 ) + f'tfU+faJ)^,*), (4) 

where cu q (t) 2 = m(t) 2 + q 2 (note that in equilibrium A^ would only depend on ti — t 2 ), and 
ti denotes the initial time of the system. The KBEs can be derived within the Schwinger- 



Keldysh formalism, see e.g. 
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5M29| ; for simplicity we assumed Gaussian initial correla- 
tions for 0. The boundary conditions for A q are fixed by microcausality and canonical 
quantization for a real scalar field, A~| tl=t2 = 0, d tl A~\ tl=t2 = — d t2 A~\ tl=t2 = 1 and 
c^c^Aql^^ = 0. The boundary conditions for A+\ tl=t2=ti are determined by the physical 
initial conditions of the system at time tj. Below, we will drop momentum indices q when 
possible. 

The quantities II 1 * 1 appearing in ()3]) and (T4]) are the self-energies of 0; in analogy to (J2]) 
they are at leading order in 0[xi] given by 



n-(xi,x 2 ) = ([o\xiM],o\xi(x2)]]) , 
n + (x 1; * 2 ) = - l -({o[ Xi {x,)io[xi{x2)\}) 



(5) 



and contain information about the interaction between <fr and the background fields Xi- They 



can be calcu 
e.g. Ref. [25|, 



ated in terms of the two-point functions of Xi within the 2PI formalism (see 



26] for details). 



III. TOWARDS BOLTZMANN EQUATIONS 

We will discuss the emergence of a description of in terms of effective kinetic equations 
by using analytical solutions of the full KBEs that are found with the WKB method. To 
this end, we make the following assumptions (we also send ti — > — oo, effects of finite ti are 
discussed below): 

1) The self-energy U ± (ti,t 2 ) is strongly damped with respect to the relative time \ti — t 2 \, 
approaching zero for \ti — 1%\ > T- mt , where we introduced the interaction time r; nt . Here, r- mt 
can be considered as definition for the duration of e.g. scattering events. It is then possible 
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to define one-sided Fourier transforms of the self-energies with respect to relative time, 

poo 

rr^w)^/ dz e iu>z I^fat-z) . (6) 
Jo 

In equilibrium the minus-component would correspond to the common retarded self-energy, 
U-(t,u) = U R (co). 

2) We assume that for fixed time t the pole structure of (oo 2 — u^(t) — n~(t,w)) _1 is 
dominated by the root oo — fl t = fi t — |T t , with 

Sk= ^ + Refi-(t,a), r t ^- lmU 2 t t ,Qt) ■ (7) 

Here, Q t and r 4 are the energy and damping rate of the ^-resonance, and we assume weak 
damping with r t <C Q t . The generalization to the case with more roots (interpreted as 
collective excitations) is straightforward. 

3) Three time scales of the system are relevant for our discussion: The above interaction 
time r int , the damping rate T t of the field 0, and the characteristic rate H ~ T t /T t ~ tl t /fl t 
with which the field changes its properties (e.g. due to a temperature change of the fields 
Xi in an expanding universe; then H denotes the Hubble rate). The main assumption, 
underlying all derivations of kinetic equations, is the separation of time scales: 

r t « rr a l and t t /T t , Q t /Q t « . (8) 

The first condition implies that the duration of one individual collision is shorter than the 
average time between different collisions; the second condition implies that the field does 
not significantly change its properties during a single collision. 

Localization of n ± (ti,t 2 ) in time occurs for different physical reasons: If all fields that 
appear inside n ± reside in a thermal bath, their propagators are exponentially damped 
with the thermal damping rate; in that case n 1 * 1 is at least damped with the same rate. 
More generally, when considering the scattering of particles, the duration of a scattering 
event typically is related to the de Broglie wavelengths of the interacting particles and is 
of the order f2 _1 ; formally the localization arises from the momentum integrals inside of 
loops. Finally, virtual particles with energies ~ M can only exist for times At < M' 1 as 
allowed by the uncertainty principle. Exceptionally small temporal damping of U^ 1 is in 
general associated with sharp features in the Fourier transform n ± (t, w) and corresponds to 
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resonant phenomena. Throughout, we will not specify the interaction time T mt , which has 
to be calculated on a case by case basis, and leave it as a free parameter. 

Under the above assumptions, the first KBE ([3]) can be approximately solved by using 



the WKB method 23 , which is formally an expansion of the solution to 03J in terms of the 
Planck constant h llj. At second order in h, we find for all times t\ and t 2 that 

sin (tdt> nAe-tift" r "l 
A-(M 2 ) = — ^ 7 =L= . (9) 

For small time separations, \t\ — t 2 \ < r int , the time-dependence of n _ (t,a;) is negligible, 
and ([9]) can be obtained by a Laplace transformation as in 15|, [32J. For \t± — t 2 \ > Tint, 
we illustrate the derivation of the WKB result by splitting A - into positive and negative 
frequency modes, A~ = f (AT^ — A7_ A Inserting the functions A^- ) (ti,t 2 ) with a = ± into 
fl3j) yields (for t l >t 2 ) 

^J(ti) - ^ - ian tl T tl ~ - / ITfa,*! - z ) e / t r^ t '(ian t ,-ir t ,) ^ (1Q) 

where we neglected terms suppressed by (jHJ) or r t /f) t ; terms containing f2 tl cancel out. Since 
n _ (ti,t') effectively vanishes for \t\ —t'\> r- mt we extended the z-integration limit formally 
to oo. Q t and T t are practically constant over the support of II - , hence the z- integral equals 
Ren(t 1 ,f2 tl ) — iaiImTl(ti,£l tl ). Then, the real and imaginary parts of ffTU]) reproduce the 
definitions of fl t and T t 0. The case t\ < t 2 is treated analogously. 

Knowledge of the spectral function A - and the self-energies n ± are enough to find a 
solution to the second KBE (jlj); it is given by the memory integral 

A + (t u t 2 ) = r dt[ f 2 dt' 2 A-fa^n+ft^A-^,*.,) , (11) 

where we allowed for a finite initial time t^. This can be obtained by using the initial 
conditions of A - and the fact that A~ solves the first KBE ([3]). 
For convenience, we now split ( TlT]) into three terms, 

A + (tx, t 2 ) = B(h, t 2 ) + C(h,t 2 ) + V(h, t 2 ) . (12) 

Here, B will correspond to the classical Boltzmann behavior, whereas C and T> will give 
corrections that can be neglected within our approximations. The terms B and C are defined 
to contain all contributions to A + that come from the time integration over t',t" < ts = 



min(ti,t2) in (II ip . where we defined the Boltzmann time t B - The term T> contains the 
remaining contributions that come from outside this region. For a self-energy that is exactly 
local in time, U + (ti,t 2 ) ~ 5(ti,t 2 ), only B and C would contribute, V would be identically 
zero. In more general cases, D is is suppressed by Y t jr^ and remains thus negligible. 

Now, B + C is split such that B contains all terms with equal-sign frequencies, B ~ 
A^A^; C contains the remaining opposite-sign terms, C ~ A7)A7_ a ). Then, C can be 
written as 

if£dt'n t , l\f t }dt'r t ,\ r t B 



/■oo 

X 



oo 







j TL+ jr, r-z ) jt_ z dt ,^ lTtl) + h n (13) 



The first term in the second line corresponds to n + (r, fi*)/f2 T . Unless this expression is 
oscillating with frequencies ±2f2 T , the integral over r is of the order Imll + (r, f2*)/f2^ and 
hence negligible. Finally, B can be written as 

cosff'^tU^e^l^^l tB 



L t 2 

coo 



X 



r dz n ^ T - z ) P s T T ^dt>^ t ,-\T t ,) , h fl4 N 



where again we obtain Il + (r, VL* T )/VL T in the second line. 
From (Tl4^) . and adopting the common definitions 

2Ren+(%) = -o t (r> + r<) , r t = r> - rf , (is) 

we obtain for the statistical propagator A + , up to terms that are suppressed by T t /fl t and 
rtA"int > 

cosff/^U'V^^'l , ts 
A+(t a ,f a ) ~ — ^ ^L= / dr(T> + T<)e-^ dt ' . (16) 



t-2 



51+2/(*b) 



We can now define the suggestive quantity /(ts) as function of t B as indicated in ({TBI) . Most 
importantly, from its very definition it follows that f(t B ) solves the BE 

d tB f(t B ) = (i + f(t B )) r< - f{t B )r> . (17) 



Then, our result for the statistical propagator becomes 

cosfr* 1 dtnA e- k *\£ dt ' r *'\ 

A + (ti,h) ~ V 2 - jj— (1 + 2f(t B )) , (18) 

2 V *1 *2 

which makes explicit that under the above assumptions the memory integral ( ITTjl is governed 
by Boltzmann behavior. Eq. (fl71) can be written as $/(t) = -T t (f(t) - /(*)), with /(t) = 
(rf /rf — If the are in local thermal equilibrium, the gain- and loss rates Tf and 

Ff are related by the detailed balance condition Tf /Tf = e~ nt ^ T( - t \ where T(t) denotes an 
effective temperature; then f(t) becomes the usual Bose-Einstein distribution. Note that 
(fT7|) and ffT8]) are valid for each field mode q separately. 



IV. DISCUSSION 



It is instructive to consider the solutions ([9]) and (|T8l) locally in time, as they appear 
e.g. in loop diagrams when calculating self-energies or when deriving properties like the 
energy density of fields. Defining mean and relative times as t = (ti + t 2 )/2 and y = t\ —t 2 , 
respectively, we obtain in the limit y H^ 1 that 

sm(yQ t ) 



A-(t!,t 2 ) 

A + (t!,t 2 ) 



-e 2 



llv|r t 



(l+2/(t))e-f 



and 



(19) 



lli/irt 



+ 2«J/(t) 



cos(y^ f 
2^ 



where = /(t) — /(^)- The damping rate of the spectral function as function of y is just 
given by IV The damping of the statistical propagator can be understood by splitting it 
up in two parts as indicated (if the background is in thermal equilibrium, they correspond 
to the equilibrium and non-equilibrium parts of the propagator, respectively). In general, 
the term proportional to Sf(t) remains undamped (cf. discussion in 30j). However, the 
overall damping of A + approaches T t if \Sf(t)\ 1. In this case, and locally in time, 
our results reproduce the common Kadanoff-Baym ansatz in the weak damping limit (see 



e.g. Refs 



.flaa 
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of a generalized phase space distribution 
Tf may include off-shell transport 



3l|). This implies that fit) as defined in (flEl) indeed plays the role 
unction of effective plasma excitations. Note that 
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33]. 



Using (TTTj) . we can discuss the effect of boundary conditions at finite time (U = for 
definiteness) . If Il + is nonsingular at t\ — i 2 = 0, ffTTl) implies the initial conditions = 
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A+ = A+ = 0; here, we defined A+ = A+| tl=ta=0 , A+ = d tl A + \ tl=t2=0 = d t2 A + \ tl=t2=0 and 
Af~ = dt 1 dt 2 A + \ tl= t 2 =o- Arbitrary initial conditions for A + can be implemented by formally 
adding SU + (t 1 ,t 2 ) = -d tl d t . 2 A + {t l ,t2)5{t l )5{t 2 ) to Ii + {t u t 2 ) in (PJ. Boltzmann behavior 
according to ( [TBI) only arises when A+ = and f^A+ = Q^Af = \ + /|t=o! otherwise, f TTTj) 
generates oscillating terms that are exponentially damped away with the rate iy 

In the early universe, quantum fields propagate in an expanding background described 
by a Friedmann- Robertson- Walker metric with scale factor a(t) and Hubble rate H = a /a. 
Since in conformal coordinates this is equivalent to a time-dependent mass term, expansion 
is already accounted for in our calculation. The BE ( TT7|) remains unchanged, but the distri- 
bution function becomes a function of comoving momentum q com = a(t)q. Going back to 
units of physical space and momentum yields the well-known result 

dtUt) = (i + / q (t)) r< - / q (f)r> + # q v q / q (t) . (20) 



V. CONCLUSIONS 

We have presented a simple method to derive BEs from the fundamental KBEs, and 
illustrated it in case of a single real scalar field <fi. Instead of performing a gradient expansion 
on the KBEs, we derived approximate WKB solutions to the full KBEs that are valid under 
the physical assumptions of weak coupling and separation of interaction and scattering time 
scales. As a result, we obtain analytical expressions for the two-point functions of — given 
by the spectral function A~(tx,t 2 ) and the statistical propagator A + (t 1 ,t 2 ) — that are valid 
at all times. Locally in time, our results reproduce the common Kadanoff-Baym ansatz. 
Although we formally neglected backreaction, it can in principle be included by reinserting 
our explicit solutions for A^ 1 into the calculation of the <ft self-energy; if the self-energies 
remain well-behaved, a self-consistent description in terms of BEs emerges. In any case, as 
long as (WKB-)solutions for the spectral functions can be found, the presented a ppr oach can 



be extended to multi-flavor problems, flavor oscillation and to fermionic fields [34J, as well 
as beyond the weak damping regime. No assumptions about the structure of distribution 
functions in Wigner-space needs to be made. Instead, the validity of a description in terms of 
distribution functions can be directly read off from a comparison between spectral functions 
and statistical propagators. 
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